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Abstract 

We study the quantum Hall states in the lowest Landau level for a single wide 
quantum well. Due to a separation of charges to opposite sides of the well, a 
single wide well can be modelled as an effective two level system. We provide 
numerical evidence of the existence of a phase transition from an incompress- 
ible to a compressible state as the electron density is increased for specific well 
width. Our numerical results show a critical electron density which depends 
on well width, beyond which a transition incompressible double layer quantum 
Hall state to a mono-layer compressible state occurs. We also calculate the 
related phase boundary corresponding to destruction of the collective mode 
energy gap. We show that the effective tunneling term and the interlayer sep- 
aration are both renormalised by the strong magnetic field. We also exploite 
the local density functional techniques in the presence of strong magnetic field 
at v = 1 to calculate renormalized As as- The numerical results shows good 
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agreement between many-body calculations and local density functional tech- 
niques in the presence of a strong magnetic field at v = 1. we also discuss 
implications of this work on the v = 1/2 incompressible state observed in 
SWQW. 
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I. INTRODUCTION 



Recent interest in the properties of high-mobility quasi three-dimensional electron sys- 
tems and double quantum well structures (DQWs) in strong magnetic fields ||f|] which 
are quite different in their structural details, has led to discovery of a novel forms of the 
fractional quantum Hall effect. In the study by Eisenstein et al. M, the sample is a tradi- 
tional symmetric double quantum well (DQWs) structure. The envelope wave function in 
the middle of the sample is small because there is a large barrier between the two wells, 
whereas in the work of Suen et al. |T|,^] the samples are actually a single wide quantum 
wells (SWQWs). The self-consistent electric field arising from the presence of the electrons 
themselves, splits the well into two spatially separated electron layers in the z-direction, 
effectively creating a double layer structure, having a small self-consistent barrier with a 
finite envelope function in the middle. 

Based on the fact that these samples are single wide quantum wells, it has been argued by 
some authors that the quantum Hall mechanism of these systems is different from DQWs. 
However, it has been claimed by Song He et al. that the result of finite size exact 
diagonalization calculation disagrees with this viewpoint and the bare Asas, the difference 
between the second and first subband energy in zero magnetic field, which is reported by Ref. 
is not the relevant parameter determining the nature of the ground state in a strongly 
correlated system. One may expect that the electronic structure of a SWQW is affected 
by a strong magnetic field and therefore it changes the Asas- In other words Asas is 
renormalized by a strong magnetic field. It has been claimed that the renormalized Ag AS is 
small enough for the Halperin state 0, to be the ground state ||. In this article we report the 
numerical calculation of the renormalized Af A5 for different realistic models of the samples 
in Ref. [[J. Our analytical expression for renormilized tunneling at v = 1 is equivalent 
with the results of Hartree-Fock approximation (HFA). We also implicate the renormilized 
tunneling's difinition for fractional filling factors e.g. v = 1/2. Our first model is a many 
body Hamiltonian calculation with electron-electron, isospin — potential interactions, where 
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one may think of the layer index as an isospin index and the second is based on the local 
density functional approximation (LDFA) in the presence of a strong magnetic field. We 
also use LDFA to evaluate the bare and renormalized tunneling at filling factor v — 1. There 
are some special cases in which SWQWs may effectively described by DQWs, i.e. high areal 
electron density and large well width. In this article DQW is the jargon for SWQW at high 
density or large width where the overlap of wave functions is negligible in the middle of 
the well. In these limits the renormalized Af AS , which is the difference energy between the 
first and second subbands, in the presence of a strong perpendicular magnetic field, is small 
enough and the quantum Hall ground states may be described by Halperin's wave functions. 
In the past few years SWQWs have been approximated by DQW, but there are significant 
correction to the DQW Hamiltonian due to overlap of electrons wavefunction which are 
localized in the SWQWs edges. In this article, we show that these corrections become more 
important either by decreasing the electron density or decreasing the well width. We also 
show that the two-component quantum Hall ground state at v — 1 which is stabilized by a 
small Af AS at high N s may evolve continuously to a one-component normal state at large 
Af A5 and low N s where we call a state of the 2DES with no quantum Hall effects by normal 
state. 

Throughout this paper we neglect any explicit consideration of electron spin, assuming 
that the magnetic field B is always high enough to totally spin polarize the system due to 
a large Zeeman splitting. We also ignore any effect of Landau level mixing and higher sub- 
bands, working exclusively in the lowest spin-split Landau level of the lowest subband. These 
assumptions are consistent with all the specific experiments to be disscused in this paper. 
It is obvious that the DQW and SWQW have some common properties which bring them 
into the same class. For the sake of simplicity, we use the term of bilayer electron systems 
(BLES) to unified DQW and SWQW as long as we concern about their common properties. 
There are two energy scales associated with a BLES in the presence of strong magnetic 
field, Ag AS and the many-body Coulomb interaction between electrons e 2 / e o^o where e is 
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dielectric constant and £o = yhc/eB is the magnetic length corresponding to cyclotron orbit 
radius with electron cyclotron energy Hu c , where uo c = eB/m*c is the cyclotron frequency. 
Competition between these two energies leads to a rich phase diagram. In a single layer two- 
dimensional electron system (2DES), the integer quantum Hall effect (IQHE) is related to 
gaps in the single particle density of states produced by the electron cyclotron energy (hco c ) , 
while in the fractional quantum Hall effect (FQHE), the gap in the excitation spectrum is 
the result of the electrons' Coulomb interaction. The quantum Hall plateau for BLES is 
controlled by the total filling factor (v tot ) which is the sum of the filling factors in each layer. 
Therefore the odd-IQHE of BLES corresponds to even denominator filling factor for each 
layer, for example v to t = 1 is produced by v = 1/2 in each layer, even though v — 1/2 in a 
single layer does not yield a plateau. Observation of odd-IQHE is due to correlation between 
electrons in different layers which is destroyed if they are uncorrelated by increasing the layer 
distance or decreasing Af AS , although for DQW with small enough separations the IQHE 
survives even in the limit Af A5 — > 0, due to a spontanous broken symmetry. This phase 
transtion has been observed for DQW by Boebinger et al. || and explained theoretically by 
MacDonald et al. @. This effect has also been observed for SWQW by Suen et al. [|IJ. In 
this article we present a theoretical study of SWQW incompressible to compressible phase 
transitions at filling factor v = 1. 

In Section |I[ we obtain the many body Hamiltonian of SWQWs in the absence of mag- 
netic field. We study the effect of strong magnetic field in Section |T| by projection of 
Hamiltonian onto the lowest Landaue level's Hilbert space and comparing with the Hartree- 



Fock approximation at v = 1. In Section [iy], we use local density functional approximation 
in the presence of perpendicular strong magnetic field to find SWQWs band structure to 
evaluate the renormalized A5^5 and comparing with other results. In Section [V], the collec- 
tive modes in quantum Hall state is studied. We present a phase diagram for quantum Hall 
compressible-incompressible phase transition and comparing with experimental results. 
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II. MANY BODY HAMILTONIAN OF SINGLE WIDE QUANTUM WELL IN 

ZERO MAGNETIC FIELD 



In a SWQW the electrons are confined in the x — y plane by an external potential barrier 
and therfore the energy spectrum becomes quantized into electric subbands. The number 
of filled subbands is a function of the areal density, N s . In experimental work at zero 
magnetic field, the first two subbands are typically filled |],@]. Neglecting high subbands, 
the subband degree of freedom SWQWs effectively reduces to a two level system. This 
system may mapped to a 2D electron system with an isospin degree of freedom. In the 
isospin language the isospin states, up (T) and down (j) refer to first and second subbands 
respectively. In the absence of perpindicular electric field bias, the system is balanced and 
the up (down) states are reflection symmetric (antisymmetric) states respectively, with bare 
eigenenergy difference A SAS . For the sake of simplicity we define As as as the eigenenergy 
differece between the two subbands in both balanced and imbalanced cases. The general 
second quantized many body hamiltonian in this picture is 

H = f + V, (1) 

where T, in analogy to the double well case is refered to as the tunneling energy 

r = -%E(W-^i)- ( 2 ) 

z p 

Here A^ A5 is the subband energy difference of noninteracting electron system in the SWQW. 
The interaction part of the many-body Hamiltonian employed in our studies can be written 
in zero magnetic field as 

v = \y,Y, vf^d&d^dpd? . (3) 

p.p'.q M 

In the above equation ty£i 'a <T1<Ta j s ^ e Fourier transform of the electron-electron Coulomb 
interaction 

y^' 2 a 1( r 2 = ( d Zl [ dz 2 i)*<(z l )^(z 2 )i) ul (z l )i) a2 (z 2 )e- q \ zl - Z2 \ , (4) 

e q J J 
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where q is the in plane wavevector. C^(CZ) is the creation (anihilation) operator for 
electrons with 2D momentum p and isospin state a and obeys, 



It is convenient to define the 4— vector isospin density operator 



(5) 



or 



(6) 



where cr° is unit matrix and a a is the ath component of Pauli matrices 
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Here /i G {0,1,2,3}. The zeroth component of is the scalar density operator and the 
rest is the vectorial isospin density. One may write the 4— vector isospin density operator in 
terms of creation-annihilation operators 



<7 r 



(8) 



The transformed hamiltonian Eq.(|I]) with continous wave vector q simplifies to 

d 2 q 



fr &SAS g , n , . 1 

H = -— a (q = 0) + 2y (2nY 



^(-q)^(q) a w (q). 



(9) 



The isospin potential V^q) in Eq.@ describes the electron-electron isospin interaction, for 
the sake of simplicity we call the V^(q), the isopotential matrix elements. They are linear 
combinations of the coulomb potential form factors in Eq.([D and are defined in Eq.([TT|). 

One finds after a little algebraic calculation, A° SAS is renormalized by the exchange 
electron-electron interaction 

d 2 q 



±SAS 



A°sas + 2 



(10) 



The enhancent may be positive or negative, depending on the sign of integral upon Vo z (q). 
For example it is positive if the density is low and negative at high density for a given well 



width. Therefore for large (small) bare A^ AS the renormalized A$as is greater (lesser) than 
A° SAS respectively. In the section |TV|, we will calculate Asas by local density functional 
approximation (LDFA) where the effect of local exchange-correlation interaction is included 
in Asas- One may show that the results of LDFA and HFA are close, hence the left hand 
side of Eq.(|H]) can be evaluated by LDFA. 

The non-zero isopotentials elements in Eq.(|9|) may be written in terms of the electron- 
electron interaction of matrix elements defined in Eq.@. 
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(11) 

Taking the advantage of symmetry properties of Eq.(f|), and the fact that the subband 
eigenfunctions are real, we define 

v n = v nn ; v ii = v iiu ; ( 12 ) 
V q n = v£ m = V^ m , (13) 

_gTJ- = yTT-14 _ priTTi _ yTUT _ y-UTT (14) 

c n = v mi = v mi = v mi = v im f (15) 

£)TI — y-UIT _ yilTJ. — y-lTU — ytlU (16) 

Note that the y component of isopotential are zero. In the next section we will see that 
these elements are still zero in the presence of the prependicular magnetic field. However, 
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it has been found that [16] the y component are not zero due to tilted magnetic field effect 



and therefore all of isopotential elements in Hamiltonian are nonzero. One may find by 
inversion symmetry of the Hamiltonian that both C^- and are zero unless the system is 
imbalanced by applying a perpendicular external electric field. 

Another extreme case may happen in the high density balanced regime, shown in Fig. ([I]). 
The overlap integral between the wavefunctions in the left and right is negligible and there- 
fore we may show that, ^(z)! = I^H- 2 )!- m this regime one finds the only nonzero matrix 
elements are B\ l , and 

yT\ = yll = yU = ytf _ ( 17 ) 

In Fig.(0), V xx (q) and V zz (q) is shown for two different areal densities. At large areal density 
V zz (q) is negligible compared with V xx (q) and the corresponding Hamiltonian is effectively 
the DQW Hamiltonian. The non-zero isopotential matrix elements in this case are Voo(q) 
and V xx (q) and may be obtained by applying Eq.(|TT|) in Eq. ([□]). These matrix elements 
have been exploited before in the study of double quantum wells [|7||| where the SWQW 
is effectively a DQW. Note that in this case Vo 2 (q) is nearly zero and the effect of many 
body exchange part of Hamiltonian on A S as is negligible so that, A S as — ^°sas- m the next 
section we will see that this is true even in the presence of strong magnetic field. Therefore in 
the DQW or high density SWQW problems, the effect of the many body part of Hamiltonian 
on A S as is zero. 

Therfore one may decompose the electron-electron interaction Hamiltonian Eq.(^) into 
three terms 

V = V DL + V B + %. (18) 

where Vdl, Vb and Vu are the double layer, balanced and unbalanced part of electron- 
electron interaction respectively 

Vdl = \Iw? Wtftf-iWq) + K(qK-(-q)^'(q)} , (w) 



Vb = ^ / {^(q)[p(-q)^(q) + ^(-q)p(q)] + K(q)fr*(-q)a*(q)} (20) 

Vu = \ I {^(q)[p(-q)^(q) + fl-q)*-(q)] + ^(q)^ (-q)**(q) + **(-q)*"(q)]} . 

(21) 

Vox has the same form as in a DQW problem. The last two terms are corrections to the 
DQW problem, due to non-zero electron density in the middle of the SWQW. These two 
terms are negligible when the density of electrons in the middle of the SWQW is small 
enough. 

III. MANY BODY HAMILTONIAN OF SWQWS IN PRESENCE OF STRONG 

PERPENDICULAR MAGNETIC FIELD 

In this section we study the SWQW in the presence of strong magnetic field where all 
electrons are located in the LLL. We exploit the techniques of projection onto the LLL which 
has been developed by Girvin and Jach |J. We will show that this technique is equivalant 
with the Hartree-Fock approximation at filling factor v = 1. 



A. Projection onto the lowest Landau level 

Taking the perpendicular component of the magnetic field to be strong, we restrict the 
Hilbert space to the lowest Landau level (LLL) and exploit the LLL projection formalism 
which was developed to study of collective excitations in the FQHE @,[l(J- The interaction 



part of the Hamiltonian which describes the low-energy excitations of the system obtained 
in Eqs.(|l8D through (|2~lD , may be projected onto the LLL by using 



P-qPq = P-qPq- Pq=O e ^ 



^qP q =^qPq-^= e^ /2 . (22) 
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In the above equations p q and are the projected total charge density and /ith component 
of isospin density operators respectively |T2| 



P(q) = 7fl e"*™ = 7fE^ ^0*) (23a) 

3 = i V ^4 j= i 



(23b) 



^(q) = ^ X>~ iq ' rj < = ^ X>~~ ^0*) 

where g = £o(q x + iq y ). The magnetic translation operator for the jth particle 

r q (j) = e^* - **** , (24) 
is a unitary operator satisfying the closed Lie algebra 

r q T k = T q+k ei qAk , (25a) 
q /\ h 

[r q , r k ] = 2i r q+k sin — — , (25b) 

where q A = £p(q x k) • z. 

One may find the commutation relations between isospin density operators by using 

[^ 1 ,^J = (e fcl ^ /2 -e fc ^/ 2 K 1+fc2 , (26a) 

K><) = ~ ek2kV2 ) ^ + ( efclfc2 * /2 + ^ K ' 2 ) n i+k2 , (26b) 

where k* is the complex conjugate of k. Note that 4-components are labeled by Greek indices 
and Latin indices denotes spatial components. The tunneling term is renormalized by the 
last linear term of Eq . (f2"2|) . The projected tunneling term is thus 

T = -^p^(q = 0), (27) 

where Ag AS is the renormalized energy difference between up and down isospin states. Due 
to the effect of the strong magnetic field on the exchange energy, we have 
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Af AS = A SAS + 2 / |^(q) e~W 2 . (28) 

One may obtain the analytical expression for renormilized tunneling, Eq.(^), using the 
Hartree-Fock calculation where the many body ground state wave function may be written 
as a single Slater determinant. A particular class of single Slater deteminant at v — 1 in the 



Landaue guage which yields to Eq . (|28|) can be written in the form [|12 



i*> = n(<Ut cos ^ir + ^ sin ^ eMX) )i°) ( 29 ) 

where |0) is fermionic vacuum and C x + . creates an electron in the Symmetric (Antiymmet- 
ric) or Left (Right) in orbit <px, respectively. Taking the many body wave function Eq.(f29|) 
and using the HFA, leads to Eq.(^). We may use Eq.(^) for any fractional filling factors 
e.g. v = 1/2 to evaluate the renormalized tunneling. Note that this generalization is just a 
simple extrapolation of Eq.(^) to fractional filling factors. As one may see from Eq . (|28|) , for 
a given well width and areal density, Vo 2 (q) is given and Af A5 is just a monotonic decreasing 
function of magnetic field 

A* AS (v) - A SAS cc TV' 112 , (30) 

where v is the Landau level filling factor. The sign is specified by the integral in Eq. (fZB]) , 
which depends on the electron areal density. Our numerical calculation shows that in the 
domain of experimental densities, the sign is minues, i.e. the renormalized As as is l ess than 
the bare one. It is instructive to write Af^ 5 (z/) in terms of the renormalized As as at given 
filling factor i.e. v = 1 and the bare Ag^s 

Af^(^) = (1 - v- 1/2 )^sas + v- 1/2 A B SAS (v = 1) • (31) 

We may use Eq.(|31|) in order to find Af^ s (z/ = 1/2). To the best of our knowledge there 
is no reliable measurement or calculation to find A B SAS and interlayer seperation (d) at 
fractional filling factor. Our numerical results for realistic samples has been shown in Table 
I. According to thses results, we may conclude that the renormalized tunneling term at 
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v = 1/2 is large enough that 331 Halperin's wave function, which is exact in the limit of 
vanishing tunneling, is not the best variational wavefunction to describe the ground state. 
However, it is not obvious that Ag AS (v 7^ 1) is the relavant parameter to the FQH gap 
corresponding to the FQHE at v ^ 1. 

At the end of this section and for further calculations we list the projected isopotentials 
onto LLL 

V = Vdl + Vb + Vu , (32) 
Vdl = \ J ^ M(q)p(-q)p(q) + V.(q)o*(-qK(q)} , (33) 

Fb 4 / + ^(-^^(^l + ^(^(-ql^fq)} > (34) 

Vu = \ I {^(q)|?(-q)^(q) + o*(-q)]?(q) + ^(q)^ (-q)^(q) +^(-qK(q)]} 

(35) 

As we mentioned in the zero magnetic field case, the last two terms of Eq. (|32|) , are negligible 
when the density of electrons in the middle of SWQWs is small enough, and the system is 
effectively a DQW with renormalized Af AS (which is nearly equal to the bare one) and 
hence one may expect that the many body quantum Hall ground state is described by the 
Halperin variational wave functions as it is in DQWs. 

IV. LOCAL DENSITY FUNCTIONAL RESULTS FOR SWQWS IN STRONG 

MAGNETIC FIELD 

The system we study numerically in this section is a wide, single, GaAs quantum well 
in the presence of a strong perpendicular magnetic field. We consider a slab of GaAs of 
thickness W confined between two infinitely high barriers, which represent Al x Gai_ x As. 
This system is doped by Si delta layers, separated symmetrically, by distance h from the 
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GaAs and which have ionized donor concentration N s per unit area. We apply the local 
density functional approximation (LDFA) to find eigenvalues and eigenfunctions of the many 
body system. By this technique one may find the charge distribution function, eigenenergies 
and the effective potential, which are calculated by self-consistent Kohn-Sham equations in a 
strong magnetic field, i.e. by solving the Poisson and Schrodinger equations simultaneously. 
In the presence of a strong magnetic field, we restrict our attention to the lowest Landau level 
(LLL), for which the single-body wavef unctions, in symmetric gauge and polar coordinates 
are 

<Pm(r,0) = 1 = , (f ) m exp (jPje im6 , (36) 
The full one-body wavefunction corresponding to jth subband may be written 

$ i>w (r, 9, z) = ip m (r, 6)^(z) , (37) 
with the corresponding charge distribution function 

nW = EElMr,f)l>iW| 2 - (38) 

j m 

In the strong magnetic field regime, where all electrons are accomodated within the 
LLL and execute cyclotron orbits with a common kinetic energy, the summation over m 
yields ttiz independent of r, the in-plain position of electrons. In the LLL where the filling 
factor is a fractional number between zero and one, only the first subband is filled, and the 



ground state can be described as an isospin ferromagnet, i.e., a phase coherent state [JTT 
In zero magnetic field under typical experimental circumstances the second subband 
is partially filled depending on areal density. In a sufficiently strong magnetic field, the 
Landau level degeneracy is high enough that electrons which were in second subband in zero 
magnetic field, will be located in the first subband. The density distribution function in 
Eq.(|38"D then reduces to 

n(z) = ^\Mz)\ 2 ■ (39) 
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Following a procedure similar to that used for GaAs/Al x Gai_ x As heterojunctions at 
zero magnetic field |TJ| , we solve for the quantized energy level of j'th subband, Ej and its 



corresponding envelope function ipj(z) satisfmg the following Schrodinger equation 

[^^ + V eS (zU(z) = E j i; j (z), (40) 

Here m* is the electron effective mass of GaAs and V e s(z) is the one electron effective 
potential of Kohn-Sham local density functional theory (LDFT) which splits into three 
different contributions 

V eS (z) = V b (z)+V h (z)+V£(z). (41) 

The quantity Vf,(z) is the built in potential due to the infinite barrier and Vh(z) is the 
Hartree term due to all coulomb interactions between electrons in the presence of a uniform 
density of background ions 

V h (z) = -^j- J dz'\z - z'\n(z') . (42) 

Vj£(z) is the exchange-correlation potential in the Kohn-Sham local density functional ap- 
proximation. The LDFA exchange-correlation energy, which has been studied by Vignale 
and Rasolt [19(1, is a functional of the scalar charge density and vector current density. In 



our approximation, we ignore the current density term in the exchange- correlation energy 
and we assume that all electrons are spin polarized along the magnetic field. The best 
functional for our purposes is the Vosko-Wilk-Nusair exchange-correlation energy, Vj](z), 



which is parametrising the Ceperly- Alder Monte Carlo calculation ||21|| , and which has been 
applied to three dimensional itinerant ferromagnets . 

Results of our numerical calculations for A$as for v = 1 filling factor and several N s are 
listed in Table I. In this Table the parameter d, which we define as the distance between 
the peaks in the charge- density profile as shown in Fig.(|5|), decreases as N s is lowered. The 
values of N s which we used have been taken from Suen et al. 's experimental data 0. At 
this point, we may compare the results of the strong magnetic field LDFA and calculation 
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of the renormalized Ag AS of the many body Hamiltonian Eq. ( f28|) , in which the tunneling 
term is enhanced by electron-electron interactions. We will evaluate the integral in Eq.(f28"D 
by knowing the Vo 2 (q) which has been defined in Eq.(|TT]). We use the zero magnetic field 
orthonormal wave functions of the first and second subband to find the Vo z (q) and hence 

^sas- 

We used the self-consistent LDFA symmetric-antisymmetric wave functions for zero 
magnetic field to calculate Af A5 numerically. The form of exchange-correlation potential 
which we chose is Hedin and Lundqvist |22[ which has been used in an investigation of 



GaAs/Al x Ga!_ x As hetero junctions in zero magnetic field [p!q] . 

One may define the relative difference between the bare and renormalized tunneling term, 
5 = AsA ^ s ^ s SAS which is a quantity to measure the deviation of SWQWs charge distribution 
due to the strong magnetic field. 5 may calculated by the integral of Vo z (q)e~ q2 ^ 2 which is 
shown for different densities in Fig.(|3]). The result of integration is positive (negative) at high 
(low) densities and depends on the long wavelength oscillations of V 0z (q). Our numerical 
calculation shows that 5 is small and even negative at low densities and becomes positive 
with increasing the density. For example, at N s = 0.5 x 10 11 cm -2 , we find 5 = —0.1 which 
is relatively small due to the long wavelength cancellation of Vo 2 (q) upon integration, and 
therefore the electron density n(z) is not affected enormously by the strong magnetic field. 
Note that the sign of S is negative, which means that the tunneling term is increased by the 
strong magnetic field. In Fig.(f|), V 0z (q) is plotted for N s = 0.5 x 10 n cm~ 2 and W = 680(A). 
In Fig.([|) we compare the electron density of a SWQW in zero and strong magnetic field by 
LDFA calculations. This comparison for N s = 0.5 x 10 n cm -2 and W = 680(A) shows that 
the two densities are very close which is indicates qualitative agreement between the LDFA 
and many body calculations. 

At very large densities Vo z (q) is a smoothly negative but close to zero and hence 5 is 
an small negative number. It is important to note that S is very small in the DQW regime 
which is equivalant to large SWQW densities and is strictly zero for the true DQW case. In 
the domain of densities which is of interest for us, both calculation of LDFA and the many 
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body Hamiltonian show a decrease of Ag^s due to the strong magnetic field on the order 
of 1 degree of Kelvin. We conclude that the effect of strong magnetic field on the electron 
density n(z) is negligible at very low and very high areal electron density. At intermediate 
densities, the effect of the strong magnetic field is significant and the electron distribution 
along the indirection is changed. One ought to be able to measure this effect experimentally. 

In Fig.(|]) we plot the result of a self-consistent LDFA calculation at strong magnetic 
field for iV s = 3.1 x 10 n cm~ 2 where the difference of charge density in zero magetic field 
and strong magnetic field is relatively large. In strong magnetic field, the charge density 
is proportional to square of the lowest subband wave function. Therefore the result of 
nonzero parameter d is due to tunneling effect through V e ^{z) barrier where the eigenenergy 
is less than the top of effective potential. Obviously in the presence of strong perpendicular 
magnetic field, d is a monotonicaly decreasing function of N s , and at a critical density N sc 
and critical width W c , the parameter d continously goes to zero. For example our numerical 
calculation shows this transition between two-layer and one-layer quantum wells occurs at 
N s = N sc = 0.5 x 10 n cm- 2 for a SWQW's width of 680(A). To the best of our knowledge, 
there is no exprimental results to report the value of A§as and d [ZS| to check with our 
numerical results. 



V. COLLECTIVE MODES EXCITATIONS OF SWQWS FRACTIONAL 

QUANTUM HALL STATE 

Recently, there has been much interest in experimental work in DQWs and SWQWs |l],f|, 
a remarkable effect is observed, namely the absence of certain IQHE at odd v for sufficiently 
small As as and large interalayer distance (d). These experimental observations have been 
explained in a phase diagram proposed by MacDonald et al. [7| for the presence or absence of 
the IQHE for different DQW parameters, associated with the loss of the isomagnon excitation 
gap. We generalize their model to SWQWs with the help of experimental parameters used by 
Suen et al. 0. In the previous sections we found that in the presence of strong prependicular 
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magnetic field only the first subband is occupied where all electrons are in the isospin 
state, up. One may expect that the low lying excitations are isospin waves (isomagnons) 
corresponding to a single flipped isospin which propagates through the system. Therefore 
we may describe the isomagnons as a bound state of an electron of one isospin with a hole 
of opposite isospin. The Hamiltonian describing the low-energy excitation of the SWQW 
system has been obtained in the preceding sections 

/\B i „ Al. 

2 ~ zv " ^ ' 2 J (27T) 1 
where the isopotentials V M!/ (q) are defined by Eq.flnp. The calculation of the isomagnon 
collective-mode energy associated with the first to second subband excitation, which is a 
type of isospin wave in isospin space, is based on the single mode approximation |7j and 
the free boson model of Holstein-Primakoff transformation which has been exploited in the 



^(q = 0) + \ J ^ ^(-q) ^(q) *„(<!) , (43) 



study of ferromagnetic spin waves problem |13| . The normalized isomagnon wave function 
is given by 

|*-(k)> = ^3-(k)|*o>. (44) 



Here I^q) = ®j | f)i is the fully isospin polarized ground state of the system in which all 
electrons are located within the first subband and S±(k) = (a x (k) ±za 2/ (k)J/2 are isospin 
lowering (raising) operators. In the absence of an external perpendicular electeric field 
where the system is balanced and the first (second) subbands are symmetric (antisymmetric) 
respectively, the full Hamiltonian can be separated into terms which conserves Sl ot {H ) 

H = T + V + H', (45) 

where T is the tunneling term and 

V ° = \ I ^o(q)p(-q)p(q) + ^(q)^(-q)^(q) + 

K,,(q)[^(-q)p(q) +p(-q)^ 2 (q) ] ) , (46) 



and a term H' which creates (annihilates) a pair of isomagnons and hence changes S 1 " 1 by 
2: 



H' = \J 7^K(q)( 5+(-q)^+(q) + 3_(-q)3_(q) ) . (47) 



2 7 (2?r) 



Using the single mode approximation |14| and commutation relation Eq.(^) yields 



(*_(k)|# |*-(k)> = £ + £(k) • (48) 

Here £ = (^o|-^o|*o) and 

e(k) = Af A5 + ^(k) e ~W _ J ^y / V te (q)( fc(q) + 1 



+ / ^ ( K(q)Mq + k) - [ V (q) + V,(q) ]h(q 

+ / ^ ( [ ^o(q) + K(q) ]/,(q) e ^ q -(-k) ) j (4g) 

where h(q) = /i(q) exp(£gg 2 /2) and /i(q) is the Fourier transform of the pair-correlation 
function. For v — 1 this is given by /i(q) = — exp(— ^q 2 /2). These results generalize 
the results which were obtained previously for the special case V z (q) = Voz(q) = for 
DQWs. The physical interpretation of Eq . fl49|) has been explained in the study of DQW 
system by MacDonald et al. 0. The first term of Eq.(^) is the dressed excitation energy 
from the first to the second subband, the second term is the Hartree correction due to the 
isospin dependent part of the electron-electron interaction and the third term is the Hartree 
enhancement of A§as associated with isospin coulomb interaction and is zero for v = 1. 
The fourth term is the difference between the self-energy of electrons in the second and first 
subbands and the last term is the energy associated with the electron-hole interaction. 

The full Hamiltonian Eq.(pE3|) may be approximated by a soluble quadratic bosonic ef- 
fective Hamiltonian by application of the Holstein-Primakoff transformation 

H b = E + e(k) S+S k + ±A(k) (6 k S_ k + b{bl k ) ) , (50) 

where b^ib^) are bosonic creation (annihilation) operators. Here H and H' are transformed 
to quadratic boson pieces Z)k £ (k) b]J>k and \ J3 k A(k) (6k^-k + fyL^-k)- The latter creates or 
destroys an isomagnon pair with opposite momentum and is responsible for broken global 
gauge symmetry. 
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One may obtain A(k) by fitting to the exact expression for the matrix element of H' 
between the zero and two isomagnon states, which leads to 

A(k) = ^K(k) e -W 2 + / VMMd + k)e^-(^) . (51) 

The effective Hamiltonian Hf, may be diagonalized by Bogoliubov transformation and 
the resulting isospin-excitation energies are given by 



E(k) = ^ 2 (k) - A'(k) • (52) 

The energy of the system is reduced by generating virtual isomagnon-pairs and therefore, 
pairing of isomagnons is prefered for the ground state. At v = 1 the dispersion relation 
Eq. ( |52"D has a gap at k = corresponding to Af As where the quantum Hall gap is specified 
by the local minimum of Eq. (|52"D . For zero tunneling, Eq.(|52|) describes a gapless goldestone 
mode corresponding to spontaneous global U(l) symmerty breaking in the ground state 
describing a superfluid associated with fluctuation in S 1 * ' P5| , in spite of violating the 



U(l) symmetry of the original Hamiltonian Eq.(pE5|) where the symmetry group of SWQWs 
Hamiltonian is Z 2 corresponding to o v — ■> — a y , while the DQWs Hamiltonian has U(l) 
symmetry group when V z (q) = Vo z (q) = 0. Note that Z 2 symmetry of Hamiltonian is also 
broken in the presence of the tilted magnetic field as we mentioned before. However at 
v 7^ 1 where the third term of Eq.(^) is nonzero but constant (depending on total filling 
factor), the collective modes are gapful even at zero tunneling term and there is no goldestone 
mode corresponding to spontaneous symmetry breaking while in DQW, spontaneous global 
U(l) symmetry breaking occurs at any filling factor u, generates the goldestone modes for 
zero tunneling. The collective modes are understood by the pole of isospin density- density 
response function within LLL which has been called isomagnons. In our approximation it 
corresponds to the local minimum of dispersion curve Fig.(|6|), which are gapful or gapless 
in special circumstances as we mentioned above. Note that this low-energy mode is the 
analogous to those in the Feynman theory of superfluid 4 He [|15||. One may associate the 



odd total filling factor quantum Hall effect with isomagnons gap where Af AS , the energy 
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difference between the subbands is the the relevant energy scale to compare with the in- 
plane coulomb energy e 2 /(eoA)) ; i- e - the competition between A$ AS and the coulomb energy 
yields the correspondig QHE. Therefore the phase boundary between incompressible and 
compressible states can be defined by the collapse of the gap which occurs at certain Af AS 
and d, while the corresponding gap at even total filling factors has been identified by the 
charge excitations (magnetorotons) which is infinitely large in our approximation Hl5| . One 
may generalize the above formalism for FQHE associated with a many body energy gap due 
to correlation between electrons in different layers which is destroyed if they are uncorrelated 
by increasing the layer distance. 

At this point we may study the effect of applying a perpendicular electric field on the dis- 
persion relation of collective modes by considering an imbalanced SWQW. As we disscused 
above, inversion symmetry is broken by perpendicular electric field and the first (second) 
subbands are no longer symmetric (antisymmetric). There is also an additional term in col- 
lective mode Hamiltonian, Eq.(f45|-flT|), due to the external electric field Hamiltonian Eq. (f2~T|) . 
We have to mention that, the corresponding additional term, Vy, changes 5**°* by 1. In this 
case the charge excitations couples with isomagnons and the ground state may be described 
by the virtual isomagnon-plasmon pairs. One may divide Vu into two parts, Vjj = Vrj + Vjj, 
where V$ (Vrj) can also be separated into two parts 




+(o) 



+ v r 



(o) 



(53a) 



u 




-+(*) 



+ Vr 



(53b) 



u 



u 



where we define V T 



u*y 




(54a) 




(54b) 



Again V# (V£) will be 



bosonized by the Holestein-Primakoff transformation 
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H° b = W ij(k) (o k S_ k + &l k 4 ) , (55a) 



2 k 



H z b = ~ £ C(k) ((n k - i)S_ k + Sl k (n k - ~) ) , (55b) 



where we have used the bosonic hard sphere condition [[[7]] for Eq. (|55b| ). Here a k (a k ) is the 
charge density excitation creation (annihilation) operator, which yield no pole in the density- 
density response function within LLL ||, e.g. magnetorotons with a huge gap equivalent 
to huj c and n k which is bosonic number operator. Therefore, Eq.(|55D may describe coupling 
between isomagnons and magnetorotons due to external electric field bias. 

Following the same calculation we have done for the ballanced SWQW Eq. (|5T|) , we may 
determine r/(k) (C(k)) from the exact result for matrix element of V$ (Vy) between the zero 
and two isomagnon-magnetoroton states 

i?(k) = / ^J^ " 7 V to (q)(l + h(q) + h(k) + e*5*(** k >&(k + q) 



Jqk-k 2 )/2 

+ S [ s(-k, -q) + s(k, q) - 2 ]) , (56) 



C(k) = / ^^It /2 ^M){\ + fe(q) + M k ) + e<^^h(k + q) 



(27T) 

+^ [s(-k,-q)+ S (k,q)-2]) , (57) 



e (gfe-fc 2 )/2 



where s(ki,k 2 ) is the Fourier transform of three-point correltion function 

s(ki,k 2 ) = l + l££ £ (^(kOS^W-kx-ka)^, (58) 



and -Bj(k) = e* k - r J = e ~ T~k(j), which is projected free electron wave function onto LLL has 
been defined by Eq. ( p^ - p^ ) and (...) is the expectation value with respect to ground state 
many body wave function, |^o)- It is possible to find analytical expression of three-point 
correltion function for v — 1 quantum Hall state where the many body ground state wave 
function may be written as a single Slater determinant [2~3|. The three-point correlation 
function in real space is 
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/ \ 1 - l z 2-zi! 2 _ l^3- z 2l 2 _ l z l~ z 3l 2 . „ £ V , . 

g[ri, r 2 , r 3 ) = 1 — e 2 — e 2 — e 2 + 2e 2 cos — -, (59) 

where = Xj + iyj is in-plane complex coordinate of particles and 

U = -{r\ + r\ + r|) + r x .r 2 + r 2 .r 3 + rx.r 3 , 

V = r 3 A r 2 + r 2 A n + ri A r 3 , (60) 

with Tj Ar k = Xjijk - yjx k . 

We can calculate the isomagnon dispersion relation Eq.fl5"2"D numerically after calculating 
realistic form-factors for SWQWs by using realistic first and second subband wave func- 
tions. We used LDFA techniques with Hedin and Lundquvist exchange-correlation energy 
functional in zero-magnetic field to obtain the subband wave functions. Therefore we may 
obtain the phase diagram for the v — 1 filling factor in a SWQW structure in the relevant 
d — A S as parameter space. 

In order to make a comparison with the results for DQWs 0, we associate the vanishing 
of the isomagnon excitation gap with the loss of incompressiblity of quantum Hall states 
to obtain a phase boundary for missing the odd integral QHE, as it has been defined for 
DQWs by MacDonald et al. 0. In Fig.(|6|) and Fig.(^) we show the results of our calculation 
at v — 1 for the isomagnon dispersion relation and the phase boundary for SWQW with 
different densities. At large d a SWQWs may be described effectively by a DQW. In this limit 
the phase boundary of the SWQW is the same as a DQW with the appropriate parameters. 
The main difference between the DQW and SWQW phase boundaries appears when d is 
small, which correspond to small width or density. In contrast to DQW, in SWQWs phase 
boundary, d goes to zero where As as is finite, A C SAS = AsAs(d = 0), due to critical charge 
density effect where the short range component of the isopin coulomb interaction is soften 
and the system evolve to one-component normal state. In Fig.fl7|), our phase boundary is 
compared with DQWs and the result of LDFA in zero magnetic field for a given well width. 
The latter is consistent with the reported phase transition As Fig.(|7D shows, the phase 
transition occurs at smaller As^s in comparison with the result of MacDonald et al. 
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for a given d. These results illustrate that our model leads to a more realistic prediction 
for phase trnsition in comparison with the previous models, however, our results can not 
capture the whole experimental datas which has been reported for phase boundary 
Note that the experimental results of Princeton group are based on the zero magnetic 
field measurements of A S as and calculated d by LDFA. 

VI. CONCLUSION 

We have presented the study of the quantum Hall effect in the single wide quantum wells. 
We showed that the effect of integral overlap in the middle of the SWQW has a significant 
effect to the quantum Hall states in these systems. 

We showed that the tunneling term (As^s) and the interlayer separation (d) are both 
renormalised due to the strong magnetic field. One may observe these effects experimentally. 
We also used and comparing different techniques, projected many body Hamiltonian onto 
lowest Landau levels, Hartree-Fock approximation and local density functional approxima- 
tion to evaluate the tunneling term and interlayer separation. Our results confirm that the 
tunneling term decreases by the strong magnetic field and hence increasing the interlayer 
separation. 

We found a phase boundary for integral quantum Hall phase transition at v — 1. We 
showed that our model cover the most part of phase boundary which has been reportaed 
experimentally || in comparison with pervious models 
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FIGURES 

FIG. 1. Numerical calculated electron wavefunction of zero magnetic field symmetric and an- 
tisymmetric subbands of a SWQW with W = 1000(1) and N s = 3.1 x 10 n cm" 2 . The envelop 
wavefunction in the middle of the sample is small due to tunneling through the self consistent 
potential barrier of the SWQW. The symmetric and antisymmetric wave functions are close to 
each other and DQW is a good approximation for this case. 

FIG. 2. V xx (q) (solid line) and V zz (q) (dash line) are shown for two different SWQW with 
W = 680(1), N s = 1.8 x 10 n cm- 2 Fig.(a) and N s = 3.1 x 10 n cm- 2 Fig.(b). The V 0z (q) is 
supressed by increasing the areal electron density. V zz in (a) is scaled by 10 and in (b) by 100 to 
show in the above figures. 

FIG. 3. Numerical results of Vo z (Q£q) vs Q^o is shown for three different SWQWs densities 
with W = 680(1). Varing areal densities N s = 1.8, 2.3 and 3.1 x 10 11 cm 2 are shown by long dashed, 
dashed and solid curves. The inset exhibits the calculated electron wave function of zero magnetic 
field symmetric and antisymetric subbands of a SWQW for N s = 3.1 x 10 n cm~ 2 . 

FIG. 4. Numerical results of Vq z (Q£o) vs QIq is shown for SWQWs electron areal density 
0.5 x 10 n cm~ 2 with W = 680(1). The integral of V 0z (Q£ )e~ Q2e o/ 2 is small and therefore the 
Asas enhancement due to strong B is small. The inset exhibits the calculated electron density 
along z— direction for two different regimes, zero magnetic field (dotted points) and strong magnetic 
field (dashed line) which demonstrates a quasi 3D electron system. These two densities profile are 
nearly the same and the effect of strong magnetic field on the charge density along the z— axis is 
negligible which is true for 3D electron system. 

FIG. 5. Electron density for self-consistent calculations in zero (solid curve) and strong mag- 
netic field (dashed curve) for N s = 3.1 x 10 n cm -2 and W = 680(1). The effect of strong B is 
decreasing (increasing) the probability of finding an electron in the edge (middle) of the SWQW 
which means that electrons are pushed to go in the middle of SWQW by the magnetic field. 
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FIG. 6. Calculated phase diagram for the v = 1 QHE in SWQW system. The lines denote the 
boundaries of phase corresponding to compressible-incompressible transition. These calculation has 
done for three differnt SWQW's width 400, 680 and 800 A. The phase boundary touch horizontal 
axis at A$ AS which corresponds to critical density N sc and critical width W c . A| A5 decreased 
with increasing of W. This is the crossover between one layer compressible and incompressible 
quantum Hall systems. The solid line is the calculated phase boundary of DQW for S— layers 
separated by a distance d which is the asymptotic limit of large SWQW's width, i.e. the behavior 
of phase boundary for the sample W = 800 is closer to the DQW. The inset exhibits the isomagnon 
dispersion relation for a SWQWs width W = 680 and the areal density N s = 1.8 x 10 11 cm -2 at the 
phase boundary where the isomagnon gap is zero. The energy of isomagnons at zero wave vector 
is specified by the Af A5 where at u = 1, the zero tunneling yields gapless dispersion relation. 

FIG. 7. The comparison between the realistic (triangles) for a SWQW with W = 680A width 
and DQW (solid line) phase diagrams are shown in this figure. The results of our zero magnetic 
local density functional results (circles) is shown. The integral quantum Hall phase transition has 
been reported in Asas /(e 2 f^o) ~ 0.05 at zero magnetic field. The results of our calcuation is 
closer to the experimental phase boundary than the results of DQW approximation. 
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TABLES 

TABLE I. Comparison of zero magnetic field As as which is obtained by Hedin and Lundqvist 
LDFA in zero magnetic field and renormalized Ag AS (u = 1), A^ s (z/ = 1) which is obtained by 
many body Hamiltonian and LDFA in the presence of strong perpendicular magnetic field and also 
renormalized Ag AS with filling factor v = 1/2, for a 680(A) wide well vs areal density N s . 

A^lQiW 2 ) A sas (K) Af as (v = l)(K) Af AS (v = l)(K) A B SAS (v = 1/2) (K) 
1.8 12.5 12.2 11.4 12.1 
2.3 9.5 8.2 8.5 7.7 
2.8 7.4 5.8 6.6 5.0 
3.1 6.4 4.7 5.7 4.0 
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